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Abstract 

We study the propagation of discrete solitons on chains of coupled optical 
waveguides where finite networks of waveguides are inserted at some points. 
By properly selecting the topology of these networks, it is possible to con- 
trol the transmission of traveling solitons: we show here that inhomogeneous 
waveguide networks may be used as filters for soliton propagation. Our re- 
sults provide a first step in the understanding of the interplay/competition 
between topology and nonlinearity for soliton dynamics in optical fibers. 
The discovery of solitons in optical fibers, three decades ago [1] , stimulated a huge amount 
of work aimed at using solitons for high speed communications [2]. Many experiments evi- 
denced the role of the Kerr nonlinearity in allowing for the propagation over long-distances 
of solitons in optical fibers [2,3]; however, this nonlinear paradigm has not yet demonstrated 
decisively its advantages over other more conventional signal propagation schemes. If one 
introduces a spatial inhomogeneity in the field equations by the linear compression mech- 
anism, one could stabilize the soliton propagation [4]. Motivated by this, one is lead to 
investigate the so-called dispersion-managed nonlinear Schrodinger equation [5] 
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where E(z,t) is the electric field, v the Kerr nonlinearity, z the propagation direction and 
t the retarded time. (5{z) is the term responsible for dispersion management. In Eq.(l) the 
term f3(z) plays the role of a spatial modulation of the kinetic term; upon discretization, this 
naturally leads to consider a network arranged in a non-translational invariant topology. 

In this paper we shall consider a discrete version of the nonlinear Schrodinger equation 
(1), describing an array of coupled optical waveguides [6]. Namely, we shall show that, 
by choosing a pertinent network of coupled waveguides, one may engineer soliton filters 
for discrete soliton propagation. We focus on waveguide networks, built by adding a finite 
discrete network of optical waveguides (the graph G°) to a chain (see Fig.l). By means of a 
a general criterion derived within a linear approximation for a relevant class of solitons [7], 
denoted as large-fast solitons, we show the possibility of controlling the soliton scattering, 
by suitably choosing the topology of G°. The analytical findings are in agreement with the 
numerical study of the full nonlinear evolution equation. 

The discrete nonlinear Schrodinger equation (DNLSE), for a general network of coupled 
optical waveguides, reads 

^ = -EA^i + A I E - I 2 E - ( 2 ) 

Here E n {z) is the electric field in the nth waveguide and A is proportional to the Kerr 
nonlinearity. The normalization is chosen to be J2n I E n (z) \ 2 — 1. In Eq.(2), (3 n j is 
proportional to the mode overlap of the electric fields of the waveguides n and j [6,9] and it 
is non-0 only if n and j are neighbours waveguides. If identical waveguides are arranged to 
form a chain, then Eq.(2) assumes the usual form idE n / dz = —/3 c (E n+ i+E n _i)+A \ E n | 2 E n 
where /3 n , n ±i = (3 C . If, on the other hand, the waveguides are arranged on the sites of a non- 
translational invariant network, a space modulation of the kinetic term occurs, even if (3 
does not depend on z. As an example, one can consider the waveguide geometry depicted 
in Fig.l. Of course, one can imagine and engineer a huge variety of network topologies; our 
aim here is to show that the network topology may be used to engineer a novel class of filters 
for soliton motion over spatially inhomogeneous networks of coupled waveguides. 
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It is well known that the DNLSE on a homogeneous chain is not integrable [6]; neverthe- 
less, soliton-like wavepackets can propagate for a long time [10], and the stability conditions 
of these soliton-like solutions may be analyzed within a standard variational approach [10,11]. 
Let us prepare, at z = 0, a Gaussian wavepacket: E n (z = 0) oc exp {— ( n ~~f°) + ik (n — Co)}- 
If A > (A < 0), one has a soliton solution only if cosk < (cos/c > 0). In the following, 
we assume A > and n/2 < k < n (positive velocities). For 7 ^> 1, the variational 
soliton-like solution should satisfy the relation: 

I COs/Cn I ,„ N 

Asoi ~ PcV^- — ■ (3) 

7o 

The long-time stability of variational solutions has been numerically checked. Hereafter, the 
term "solitons" denotes such variational solutions. Interestingly, discrete solitons have been 
experimentally observed in optical waveguides [12]. 

Let us consider the inhomogeneous network obtained by attaching a finite graph G° to a 

single site of a chain (see Fig.l). We denote the generic waveguide with latin indices 

The index i can be an integer number or a greek letter a, (3, . . . according to if the waveguide 
% belongs to the chain or to G° respectively. A single link connects the waveguides of the 
chain with the waveguide a of G°. We suppose the waveguides of the infinite chain to be 
identical, so that their coupling terms are set to the constant value f3 c . Soliton propagation 
in a generic inhomogeneous network is conveniently described within graph theory [13]. The 
adjacency matrix A of G° is defined as: A° ■ = f3 iy j when % and j are nearest-neighbours sites 
of G°, and otherwise. G r denotes the graph obtained by cutting the site a from G°, and 
A r is its adjacency matrix (see Fig 1). The energy levels of G° and G r are the eigenvalues 
of A and A r respectively. 

The scattering of a soliton through this topological perturbation has been numerically 
studied in the following way. At z — 0, we prepare a Gaussian soliton, far left from (i.e., 
Co < 0), moving towards n = (sin(/c ) > 0), and with a width 70 related to the nonlinear 
coefficient A through Eq.(3). We numerically evaluate the nonlinear evolution of the electric 
field from Eq.(2). The group velocity is v — (3 c smko and, when z s |£o|/sin(fco), the 
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soliton scatters through the finite graph G°. At a position z well after the soliton scattering 
(z ^> z s ), we evaluate the reflection and transmission coefficients 1Z = J2 n <o I E n (z) | 2 and 
T = Zn>o\E n (z) | 2 . 

The interaction between the soliton and the defect is characterized by two length scales 
[14]: the soliton-defect interaction length z int = •y //3 c smk and the soliton dispersion length 
Zdisp = 7o/(4/9 c sin (I/270) cos /c). When 70 ^> 1 and z int <C Zdi sp , the soliton is a large-fast 
soliton: i.e, during the scattering it may be regarded as a set of non interacting plane waves 
and the transmission coefficients may be computed by considering the transport of a plane 
wave across the topological defect. The results obtained within the linear approximation 
are in good agreement with the numerical solution of Eq.(2) [see Figs. 3-4]. We point out 
that, even if a linear approximation is used [7], the nonlinearity still plays a role: it keeps 
together the soliton during its propagation (see Fig. 2). 

In the large-fast soliton regime, the momenta for perfect reflection and transmission are 
completely determined by the spectral properties of the graph G° [7]. The linear eigenvalue 
equation to investigate is — ^2 m A n>m E m = fJ,E n , where = (5^ is the generalized ad- 
jacency matrix of the whole network. The momenta k corresponding to perfect reflection 
(TZ(k) = 1) and to perfect transmission (T(k) = 1) of a plane wave are determined by 
imposing the continuity at sites and a. One obtains 1Z — 1 if 2/3 c cos/c coincides with 
an energy level of G°, while T = 1 if 2f3 c cosk is an energy level of the reduced graph G r 
[7]. This argument can be easily extended to the situation where p identical graphs G° are 
attached to n — 0. In the limit of an infinite inserted chain, the soliton propagates in a star 
graph, which has been recently investigated in the context of two-dimensional networks of 
nonlinear waveguide arrays [8]. 

The general analysis carried out in [7] allows for an easy identification of the graph G° 
selecting the reflection (or the transmission) of a particular k (a filter). For instance, a 
transmission filter may be obtained by inserting a finite chain of 3 sites (a, f3 and 7; see 
Fig. 3). The momentum of perfect transmission k^ is determined by 2/3 c cosA;^ = —ftp,*/- 
As previously discussed, the value of the perfect transmitted momentum does not change if 
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p identical chains are attached in n = 0, but in this case the minimum becomes sharper. In 
this simple case an analytical calculation, in the linear regime, of lZ(k) gives 



A reflection filter can be obtained by adding a finite chain of 2 sites (a and (3). A 
diagonalization of A° proves that a total reflection is obtained at the momentum kS max ^> 



setting /3p n = and p — 1. In Fig. 3, the good agreement between numerical findings, 
obtained considering the full non-linear evolution, and analytical results is evidenced for 
large-fast solitons (70 = 40). 

A high-pass/low-past filter, allowing the transmission of high/low velocity solitons (see 
Fig. 4), can also be obtained. The graph realizing a high-pass filter is given by p finite chains 
of length 2 and the coupling terms fixed to f3 c . The analytical expression for 1Z(k) is obtained 
by setting f3 a> p = j3 c and /3p a = in Eq.(4) (the cutoff on the momentum depends on p). 
The low-pass effect is obtained with a linear chain of 3 sites (a, (3, and 7). If (5p tl = 2/3 c , 
then T = 0forA; = 7r/2 and lim^^Tl = 0. Therefore, one has a low-pass with the cutoff 
momenta depending on f3 a> p (in Fig.4 /3 0iQ = (5 a ^ = f3 c and (3p tJ = 2f3 c ). 

In conclusion, we showed that, by a pertinent choice of the topology of the graph G°, one 
is able to control the reflection and the transmission of traveling kinks. The results obtained 
show the remarkable influence of topology on nonlinear dynamics, and apply in general to 
soliton propagation in discrete networks whose shape (i.e., topology) is controllable. As 
these results suggest, we feel that it is now both timely and highly desirable to develop the 
investigation of nonlinear models on general inhomogeneous networks of coupled networks, 
since one should expect new and interesting phenomena arising from the interplay between 
nonlinearity and topology. 

Acknowledgments: We thank M.J. Ablowitz for stimulating discussions. 





satisfying 2f3 c cos k! <max) 



Pa ; p. The analytical expression for TZ(k) is obtained from Eq.(4) 
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FIG. 1. The inhomogeneous networks of waveguides, obtained attaching to a chain the graph 
G°, which in this case is a finite chain of length 3. In the inset we plot the corresponding graph. 
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FIG. 2. A soliton [ko = 1.8, 70 = 40 and A given by Eq.(3)] scatters through the graph of 
two sites described in Fig 3. The predicted reflection coefficient is close to 1. Inset: wavepacket 
evolution for initial momentum ko = 0.2 (wrong non-nonlinearity) , the wavepacket spreads before 
hitting G°. 
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FIG. 3. The reflection coefficient vs. the momentum k for a transmission filter and a reflection 
filter. For the transmission filter /3o, a = P a ,p = Pc and = /3 c /2, for the reflection filter /3o, a = /3 C 
and /3 a;; 3 = (3 c /2. Stars and circles are numerical results. The reflected and the transmitted 
momentum is w 1.8. Solid lines correspond to Eq. (4). 




FIG. 4. The reflection coefficient vs. the momentum k for a low-pass filter and a high-pass 
filter (and p = 7). Stars and circles are numerical results of Eq.(2). Solid lines corresponds to 
Eq.(4). 



8 



